Hopping models and ac universality 
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Some general relations for hopping models are established. We proceed to discuss the universality 
of the ac conductivity which arises in the extreme disorder limit of the random barrier model. It is 
shown that the relevant dimension entering into the diffusion cluster approximation (DCA) is the 
harmonic (fracton) dimension of the diffusion cluster. The temperature scaling of the dimensionless 
frequency entering into the DCA is discussed. Finally, some open questions about ac universality 
are mentioned. 
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INTRODUCTION 

Hopping models have been studied for many years 
in different contexts |lj ||, [|. Our main concern 
here is the modeling of ac conduction in disordered 
solids, a field where hopping models are quite successful. 
There is a wealth of ac data for various disordered solids 
like ionic conductive glasses, amorphous semiconductors, 
polymers, etc. These solids behave similarly; from bulk 
ac data alone it is not possible to distinguish ionic and 
electronic conduction. The "universal" ac features || 
are: At low frequencies the conductivity is constant, at 
high frequency it follows an approximate power law with 
an exponent which is less than one but slowly approaches 
one as frequency increases. As temperature is lowered the 
dc conductivity goes to zero and the frequency exponent 
goes to one (at any fixed frequency much below phonon 
frequencies) . The frequency marking onset of ac conduc- 
tion is proportional to the dc conductivity, while both 
are strongly temperature dependent (usually Arrhcnius). 

In this paper we first establish hopping models as part 
of a more general framework. We then proceed to dis- 
cuss the extreme disorder limit of hopping models. In 
this limit the ac conductivity in properly scaled units 
becomes independent of the jump frequency probability 
distribution. Finally, new results are established for the 
diffusion cluster approximation, a recently proposed an- 
alytical approximation for the universal ac conductivity. 



DIFFUSION AND HOPPING MODELS 

First, we briefly consider diffusion in general and define 
the wave-number and frequency-dependent diffusion con- 
stant. Then we specialize to hopping of non- interacting 
particles on a cubic lattice and derive a number of re- 
lations. Many of the points made below can be found 
in the literature, but to our knowledge they have never 
been concisely written down aimed at hopping models. 

If p is the particle density and J the particle cur- 
rent density, the most general diffusion-type constitutive 
equation obeying linearity, causality and space-time ho- 



mogeneity is: 
J(r,t) = - / dr 



dt' D(r',t')Vp(r-r',t-t'). (1) 



It is convenient to introduce the so-called Laplace fre- 
quency s given by 

s = iu> , (2) 

where u> is the ordinary frequency. The wave-number 
and frequency-dependent diffusion constant, D(\t,s), is 
by definition the Fourier-Laplace transform of D(r' ,t'): 



dr' e~ lk r ' 



dt' e- st ' 



D(r',t') 



= I dr' e- lk r 'D(r',s). 



(3) 



It is possible to express Z?(k, s) in terms of equilibrium 
density fluctuations. To show this we first Laplace trans- 
form the equation of continuity, p = — V • J, and calcu- 
late the right hand side by means of Eq. ([!]) . If p is the 
Laplace transform of p this leads to 

sp(r, s) - p(r, 0) = J dr'D{r',s)\7 2 p(r -r',s) . (4) 

We proceed by Fourier transforming this equation. If 
Pk is the Fourier transform of the density, pk = 
J drp(r) exp(— ik • r), one finds 



spk(s) - /9k(0) 



-D(k,s)k 2 p k ( S ) 



(5) 



Equation implies pu(s) — Pk(0)/ [s + D(k, s)k 2 ] . 
Multiplying this relation by p_k(0) we find after ther- 
mally averaging 



(Pk(a)p- k (0)) = 1 
(Pk(0)p-k(0)) s + D{k,s)ki 



(6) 



We are mainly interested in the k — > limit. The 
frequency-dependent diffusion constant D(s) is defined 
by 



D(s) = lim DCk.s) 



(7) 
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We proceed to derive an equation relating D(s) to the 
mean-square displacement as function of time, consider- 
ing from now on only the case of non-interacting parti- 
cles. If d is space dimension and Ar(£) is particle dis- 
placement in time t, the relevant expression is (with an 
implicit convergence factor lim e ^o exp(— et), e > 0): 



2 roo 

D(s) = Yd], < A r 2 W>e- st *. 



(8) 



The proof proceeds as follows. If the par- 
ticle positions are denoted by rW we have 
/°k(i) — Ylj ex P [ — ^ ' r ^ CO] ■ Non-interacting par- 
ticles move uncorrelated and thus ((0k(OP-k(O)} = 
^■(exp [— ik ■ ArW (£)])■ Consequently, it is enough to 
consider the motion of just one particle and 



(PkWP-k(O)) 
(pk(O)p-k(O)) 



-ik-Ar(t) 



(9) 



Being only interested in the limit of small k, the exponen- 
tial on the right hand side is expanded: exp (— ik • Ar) = 
1 — ik • Ar — |(k • Ar) 2 + .... For the thermal average we 
have by symmetry (Ar) = and thus (exp (— ik • Ar)) = 
1 — g((k • Ar) 2 ) + .... Spherical symmetry implies that, 
if Ari is the i'th coordinate of Ar etc, (AriArj) = 
{Ar^Sij/d. Thus ((k • Ar) 2 } = (hAnkjArj) = 



k 2 {Av 2 )/d. Substituting this into Eq. (| 
forming, and utilizing Eq. (^|) leads to 



Laplace trans- 



1 



1 



L>(k,s)fc 2 



(10) 



For small k the left hand side is expanded as follows: 
s _1 (l — s _1 I?(k, s)k 2 + ...). Comparing this to the right 
hand side and letting k — * proves Eq. da). 

The case of "ordinary" diffusion is characterized by 
(Ar 2 (t)) = 2dDt where D is the ordinary diffusion con- 
stant. Substituting this into Eq. (||) gives D(s) = D at 
all frequencies. Thus the definition of D(s) is a consistent 
generalization of the ordinary diffusion constant. Finally, 
we note that Eq. (^) may be rewritten by performing two 
partial integrations: 



1 [°° 



{Q)-v(t))e- st dt. 



(11) 



Henceforth we specialize to the case of non-interacting 
particles hopping on a finite cubic lattice in d dimensions. 
Hopping models have been discussed for several years in 
the literature. Some reviews of hopping are listed in Rcfs. 
0, |[ ||, 0], some recent interesting papers dealing with 
ac hopping are given in Refs. §, 0, |]. In order to be 
able to speak about well-defined displacements we shall 
not impose periodic boundary conditions. This implies 
D(0) = 0. To calculate D(0) in the physically realis- 
tic "bulk" limit one should calculate D(s) at a non-zero 
value of s, let lattice size go to infinity and only there- 
after let s go to zero. Alternatively, one imposes periodic 



boundary conditions and doesn't consider positions at all 
but velocities, and calculates D(s) by means of Eq. jll| ) 

!■ 

As mentioned already, it is only necessary to consider 
the motion of one particle. A hopping model is defined 
as follows. Adopting the bra-ket notation of quantum 
mechanics the lattice sites are denoted by Any state 
of the system \P) is given by an expression of the form 
\P) = J2 m Pm\ m )> wnere Pm is the probability to find 
the particle at site m. The state space is a real Hilbert 
space when equipped with ordinary inner product but, of 
course, only states with positive probabilities summing 
to 1 are physical. Time development is described by a 
master equation which is conveniently written in terms 
of a "Hamiltonian" as follows: 



d 
~dt 



H IP) 



(12) 



Here, H is the well-known master equation time- 
development matrix constructed from the transition 
probabilities. H has the following properties All 
diagonal elements of H are non-positive (giving the de- 
cay rate from a state), all off-diagonal elements are non- 
negative (giving the transition rate from one state to an- 
other). The condition J^. Hj m = ensures probability 

conservation. Finally, if Pm ' is the canonical probability 
to find a particle at site m in thermal equilibrium, for 
any states \j) and \m) the principle of detailed balance 
demands that Hj m Pjh = H m jPj. 

Since the solution of Eq. (|TJ) is \P(t)) = 
cxp(Ht)\P(0)) the mean-square displacement in thermal 
equilibrium is given by 



(Ar 2 (t)> = £(r 



(j\e m \m)pW. (13) 



As always when dealing with a master equation it is 
convenient to switch to the "symmetric" representa- 
tion |^0| : First we define the operator T by T\m) = 

J P^J | to) = | Tpm). It is easy to show that the 
principle of detailed balance implies that the operator 
H = T~ 1 HT is Hermetian (real and symmetric). Fur- 
thermore, (j\exp(Ht)\m) — (j\ exp ( THT~ 1 t ) \m) = 



OlTexpCfftjT-V) = (iP 3 \eMHt)\4> m )/P 
tuting this into Eq. (|13|) leads to 



(0) 



(Ar 2 (£)) = ]T(r 3 



\>J'n 



Substi 



(14) 



From this important conclusions are arrived at. First, 
note that if any state may be reached from any other 
- "ergodicity" - H has a unique "ground state" with 
eigenvalue jl0|]. This state, |0), is the symmetric 
representation of the state of thermal equilibrium. It 



3 



is straightforward to show that this normalized eigen- 
state is given by |0) = J2 m \ip m ). All other eigenval- 
ues of H, denoted by — j n , are strictly negative (i.e., 
7„ > 0). If the corresponding eigenstates of H are de- 
noted |n), the usual trick of "sandwiching in" the or- 
thonormal set of eigenstates leads to \ exp(Ht)\ip m ) = 

Pf ] Pm +J2,Mj\ n )( n \' l l J rn) exp(-j n t). Substituting this 
into Eq. (Ful) leads to 



(Ar 2 W > = Y, " P^P^ ~ E 



jm 



■tnt 



(15) 

where \i n = -J2j m ( r i ~ r m) 2 (il>j\n}{n\il>m}- Since 
(Ar 2 (t)) = for t = Eq. ([l5|) may be rewritten as 



{Ar 2 (t)) = 5> n (l-e" 



(16) 



We proceed to prove that /i„ > 0. Suppose 
that Aj are real numbers summing to zero. 
Then J2j m ( r i _ r ™) 2 ^™ ^ °> because the 
sum may be rewritten as (j^^ Ayr^ (E m Am) + 



2 (E, 



and the first two 



(Em A m r; 

terms vanish. It follows that \x n > because for any \n) 
one has Ej^jl™) = (0| n ) = 0- 

We are now able to prove a result for the mathematical 
structure of the mean-square displacement in hopping 
models: By repeated differentiations Eq. ( |l6| ) implies 

!<Ar 2 (i)> > 



dt 2 



(Ar 2 (t)} < 
... etc. ... 



(17) 



This is the first important consequence of Eq. (|l4| ). A 
similar result applies for D(s). Substituting Eq. ( |16[ ) into 
Eq. (|) leads to 



s 2 x 



2d^ 

n 



i 

s 

SJn 



'In 



s + Jr. 



(18) 



Obviously, D(0) = and D(s) > for s > 0. Equation 
( |l8| ) is further simplified by writing in the numerator s = 
s + J n - In, leading to 



D{s) = D(oo) - i^ Mn 



In 



' S + J n 



(19) 



where D(oo) = En Mn7n/2d . It is straightforward to 
show that D(oo) = J2j m ( r i _ r rn ) 2 T jm P^ /2d where 
Tj m is the transition rate from site \m) to site \ j) (Tj m = 



—Hj m for j ^ m). Differentiating Eq. ( |l9| ) with respect 
to s repeatedly for s > leads to 



ds 



2 D(s) < 
... etc. ... 



(20) 



Substituting s = iuj into Eq. (|19|) immediately gives a 
result first derived by Kimball and Adams [jllj : The real 
part of the diffusion constant is always an increasing func- 
tion of u>. 

Our final general result relates to the case s > 0, where 
the approximate power law exponent obeys: 



< dln_D < . 
— dins — 



(21) 



The first inequality simply expresses the fact that -D(s) 
is an increasing function of s (Eq. (pp|)). Now, if 
fi(s) and /2(s) are two positive functions each obey- 
ing din //dins < 1 then this condition is also obeyed 
for fx + f 2 : We have sf'Jfx < 1 or f[ < fx/s and 
similarly for f 2 . Adding these two inequalities leads 
to f[ + f 2 < (fx + h)/s or dln(/i + h)/d\ns < 1. 
Next, we note that the function /(s) = Cs/(s + j) obeys 
< din //dins = j/(s + j) < 1 whenever C > and 
7 > 0. The second inequality Eq. ( pl| ) now follows from 

Eq. m. 



UNIVERSALITY OF THE AC CONDUCTIVITY 
IN THE EXTREME DISORDER LIMIT 

If the particles are charged the fluctuation-dissipation 
theorem tells us that the frequency-dependent electrical 
conductivity <r(s) is proportional to D(s). In particular, 
in hopping models there are relations for <j(s) analogous 
to Eqs. (|20| ) and (pl|), and the real part of the conduc- 
tivity, a', is an increasing function of to. 

We now specialize to hopping models with only 
nearest-neighbor jumps, and focus on the simplest model 
of this kind, the random barrier model (also called the 
symmetric hopping model). For this model all sites on 
the cubic lattice have equal energy. Thus, detailed bal- 
ance implies for the jump rates Tj m = T m j. Because 
we want to model a disordered solid the jump rates are 
taken to vary randomly and uncorrelated from link to 
link. Moreover, we shall assume that the jump rates are 
given by a free energy barrier E: T — T exp(—(3E). Here 
(3 is the inverse temperature if classical barrier hopping 
is considered and inverse wave function size if quantum 
mechanical tunneling is considered. The model is com- 
pletely defined in terms of the free energy barrier proba- 
bility distribution, p(E). 

At large values of (3 the jump rates vary many orders of 
magnitude. The limit p — > 00 is termed the extreme dis- 
order limit. In this limit, although all jump frequencies 
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go to zero and consequently a — > 0, one may ask how the 
conductivity relative to the dc level, a = £r/cr(0), behaves 
as a function of frequency [l3| [l4| . It turns out that a 
as function of a suitably scaled frequency in the extreme 
disorder limit becomes independent of the barrier dis- 
tribution. This result, "ac universality," applies for any 
smooth p(E) which is non-zero at the percolation thresh- 
old. More accurately, ac universality means that for any 
fixed scaled [Laplace] frequency s, as (3 — > oo the conduc- 
tivity a(s) converges to some value which is independent 
of p(E). In computer simulations we find that the larger 
s is, the larger j3 must be before there is convergence. 
Since s ~ 1 defines the region below which the conduc- 
tivity is virtually frequency independent, this means that 
the convergence to universality is fastest around the on- 
set of ac conduction. 

Although no mathematically rigorous proof of ac uni- 
versality exists, there is convincing evidence for ac uni- 
versality from several sources ||. First of all, ac uni- 
versality is clearly seen in computer simulations, where 
several quite different p(E)'s lead to the same a(s) in the 
extreme disorder limit. Secondly, the effective medium 
approximation (EMA) implies ac universality in the ex- 
treme disorder limit with the following prediction fl4| : 

a In 5" = s. (22) 

This equation was first derived by Bryksin as the EMA 
solution of the hopping model describing electrons tun- 
neling between random positions in space fll5[ . Finally, 
it is possible to physically understand the origin of uni- 
versality; basically it stems from the fact that percolation 
dominates conduction in the extreme disorder limit 

The EMA universality equation (^2j) has the correct 
qualitative features of the universal ac conductivity: a 
constant low frequency conductivity and an ac conduc- 
tivity which at high frequencies follows an approximate 
power law; the exponent is below one but converges loga- 
rithmically to one as frequency diverges. Quantitatively, 
however, Eq. ( p2|) is not accurate in three dimensions, 
and even less accurate in two dimensions. In both cases 
the onset of ac conduction is less dramatic than predicted 
by Eq. (|22|). One may speculate, however, that Eq. j2^ ) 
becomes exact in and above 6 dimensions, because here 
mean-field theory for percolation is exact (when regarded 
as a critical phenomenon), and EMA is a sort of mean- 
field theory. We have preliminary simulation data for 4 



dimensions indicating that Eq. (22) indeed works better 
as dimension is increased. 

Before proceeding to a discuss a more accurate analyti- 
cal approximation to the universal ac conductivity, let us 
briefly sketch how Eq. (|22|) is derived because this is rel- 
evant for the following. The units used are "rationalized 
units" where conductivity, diffusion constant and jump 
rate are all identical in the ordered (homogeneous) case. 
The general EMA equation for finite disorder is a rather 
complicated self-consistency equation for a(s) 0. This 



equation involves the diagonal element of the Green's 
function for a "homogeneous" random walk with uniform 
jump frequency T, G = (r|Go|r), where Go = l/(s — H ) 
is the resolvent operator for the uniform-case "Hamilto- 
nian" Ho. Gq is a function of T(= cr), besides of course 
also a function of s. In the extreme disorder limit the 
EMA self-consistency equation reduces [H to 



In cr oc sGq . 



(23) 



It is straightforward to show that sGq oc s/a for small 
s in two and more dimensions p^ j. When this is substi- 
tuted into Eq. (^3|) one arrives at Eq. (^2|) after a suitable 
rescaling of s to define the dimensionless quantity s: 



/(/?) 



a(0) 



(24) 



where /(/3) oc (3. 

As mentioned already Eq. ( |22| ) does not give an ac- 
curate representation of the universal ac conductivity in 
three dimensions. A better fit to data is provided by what 
we term the "diffusion cluster approximation" (DCA), 
which leads fllql to 



In a 



g\ do/2 



(25) 



Our simulations in 3-d are well fitted by do = 1.35. An 
equation similar to Eq. (|2^) was derived by Zvyagin in 
1980 17 by reference to cluster size statistics at perco- 
lation. Zvyagin's equation, however, is real: He has u> 
where we have iuj and Re(<r) where we have a. 

In the original derivation of Eq. ( p5|) dg was identified 
with the dimension of the "diffusion cluster" fl6| . This 
set is defined as follows. We first recall that at extreme 
disorder conduction takes place on the percolation clus- 
ter 0, [H| HFJ. More precisely, links with jump rates 
much smaller than the "percolation jump rate" can be 
removed without affecting the overall (dc) conductivity 
of the lattice; for any finite (3 this leaves us with the "fat" 
percolation cluster 0] (which becomes the true percola- 
tion cluster as /3 — > oo). However, more links may be 
removed. First of all, those on dead ends of the percola- 
tion cluster may be removed, leaving the so-called back- 
bone. Moreover, if there are two different paths between 
any two points below the correlation length, because of 
the extreme disorder one of them is much more favorable 
than the other which may be removed. After this diluting 
of the percolation cluster the remainder is by definition 
the diffusion cluster. Now, our basic assumption is that 

• In the extreme disorder limit not only dc but 
also ac conduction takes place mainly on the 
diffusion cluster. 

This assumption runs contrary to the traditional use 
of percolation theory to calculate cr{s). The assump- 
tion should be checked by computer simulations, but we 
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have not yet done so. Please note that, at any finite (3 
the assumption applies only at not too high frequencies; 
there are always isolated islands of well-conducting re- 
gions which are important at sufficiently high frequencies. 
Our conjecture, however, is that the scaled frequency be- 
low which the conjecture applies diverges as (3 — > oo. 

What is the structure of the diffusion cluster? Let 
us first recall the "nodes-links-blobs" model of the back- 
bone. According to this model [pi] the backbone com- 
prises links (i.e., quasi one-dimensional strings) and 
nodes at the intersection of links. The "blobs" refer to 
the fact that for — 1 percolation there are occasional 
"strongly bonded" regions along any link. However, as 
argued above, blobs are unimportant in the extreme dis- 
order limit of a continuum distribution of jump rates be- 
cause, of any two different paths between two points one 
will strongly dominate. 

Our picture is now the following: At any finite /?, 
the time scale on which the particle moves more than 
the node-node distance corresponds to frequencies where 
the conductivity is frequency independent; here we agree 
with the standard use of percolation theory to ac phe- 
nomena. On a smaller distance scale the diffusion cluster 
is fractal. In our original derivation of Eq. (25) we used 
the EMA equation Eq. (^3|) for hopping on the diffusion 
cluster and identified e?o with the dimension of the dif- 
fusion cluster, having in mind the fractal dimension d. 
More correctly, d should be identified with the so-called 
harmonic or fracton dimension d, which for homogeneous 
(uniform jump rate) random walks on the diffusion clus- 
ter gives the probability P{t) to be at the same place as at 

t = 0: By definition P(t) oc t~ d / 2 and Go is basically the 
Laplace transform of P(t), which implies sGq oc (s/a) d ^ 2 . 
The connection between d and d [Ell, E2| is 



d = 



2d 
2+7' 



(26) 



where the exponent 5 is given by (Ar 2 (t)) oc i 2 /( 2+l5 ) for 
a homogeneous random walk on the diffusion cluster. 

Is it possible that d = d, corresponding to 8 = 0? In 
our opinion the answer is yes, because 5 = for any 
fractal without dead ends and without loops, e.g., a self- 
similar curve like the Koch curve 0, |2^| . If the links of 
the nodes-links model for the diffusion cluster are of this 
kind one has d = d. Indeed, we argued above that the 
diffusion cluster has no dead ends or loops. 

We finally briefly discuss the temperature scaling of 
the dimensionless frequency in the DCA. As mentioned 
above, for homogeneous random walks on the diffusion 
cluster sGq oc {s/a) d °/ 2 where d a = d. When this is 
substituted into Eq. ( p3| ) we find that the function /(/3) 
in Eq. (g) is given by f{fi) oc f3 2 / d °. For d = 1.35 one 
gets / oc /3 1 ' 48 . In our simulations we find a similar power 
law for f{0)i albeit with exponent 1.37 0]. Roling finds 



an exponent equal to 1.3 



OPEN QUESTIONS 

Many questions about ac universality in the extreme 
disorder limit remain unanswered: 

• For the random barrier model: How does the uni- 
versal ac conductivity depend on dimensionality? 
In particular: Is Eq. ( p2| ) exact in 6 dimensions 
and above? 

• For the more general asymmetric hopping model, 
i.e., with lattice sites of differing energies: Does 
this model also have ac universality in the extreme 
disorder limit? If yes: Does the universal ac con- 
ductivity depend on the choice of transition rates? 
How does it depend on dimensionality? 

• More realistically one should consider hopping with 
only room for one particle at each site. For "Fermi 
hopping:" Is there ac universality in the extreme 
disorder limit? If yes: Is this the same as the ac 
universality for the random barrier model, which is 
the mean-field (Hartree) limit of the Fermi model 

©0? 

In our opinion much remains to be done in this chal- 
lenging field of research. Compared to the 1970's and 
1980's one now has the possibility of extensive com- 
puter simulations at hand. This is likely to bring further 
progress in the field. 
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